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1. Introduction 

It was discovered in the early seventies that certain well-known partial differential equations (e.g., 
the Korteweg-de Vries equation) have a very natural interpretation as (infinite-dimensional) completely 
integrable systems, and that the results from the classical study of finite-dimensional integrable systems 
could be used to clarify and to study the behaviour of these systems. The discovery of this relation 
also caused a revival of the study of (finite-dimensional) completely integrable systems, which has been 
undergoing, during the last decades, a number of significant developments. On the one hand, many 
new examples of (finite-dimensional) completely integrable systems were constructed (e.g., reductions 
of infinite-dimensional systems), and a lot of different techniques (e.g., algebraic geometric methods) 
were developed to study and solve these systems. On the other hand, the study of infinite-dimensional 
integrable systems revealed the existence of a number of interesting concepts (e.g., recursion operators) 
related to the integrability of a system, leading to a more "structural" study of integrable systems, 
which aims at a better understanding of integrability and of the structure of integrable systems. 

One of the concepts that evolved from this second line of research is that of a bi-Hamiltonian 
structure for a mechanical system. Roughly speaking, a dynamical system is said to be bi-Hamiltonian 
if it can be written in Hamiltonian form with respect to two different Poisson structures {•, -}i and 
{•, -}2, i-e., 

x = {x,H 1 } 1 = {x,H 2 } 2 

for some functions H\,H2 € C°°(M), under the additional assumption that the brackets and 
{•, -}2 are compatible, i.e., their sum is again a Poisson bracket. The existence of different (compatible) 
Poisson brackets generating the dynamics of a system was first observed by F. Magri [Ma] in the case 
of the Korteweg-de Vries equation, and it was remarked by the same author that the existence of a 
bi-Hamiltonian structure for a mechanical system is closely related to the integrability of the system. 
Indeed, if we suppose the Poisson bracket {-,-}i to be non-degenerate, we can define a recursion 
operator 

N = QP~\ 

P and Q being the Poisson tensors associated to the brackets {■, -}i and {■, -} 2 respectively. This 
recursion operator N can be used (as is shown, for example, in [MM]) to construct a family of 
functions F i: which are first integrals of the Hamiltonian system. Moreover, these functions turn out 
to be in involution with respect to both Poisson brackets, 

{F i ,F j } 1 = {F i ,F j } 2 = 

for all which, in the case of finite-dimensional systems, suffices to show the integrability of the 
Hamiltonian system. 
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We remark however that the precise relation between bi-Hamiltonian structures and integrability 
is not completely understood at this moment. On the one hand, starting from an infinite-dimensional 
bi-Hamiltonian system, the recursion operator construction does not always yield a sufficient number 
of first integrals to prove complete integrability. Also, in the case of degenerate systems, (i.e., systems 
whose associated Poisson structure is degenerate,) the recursion operator does not exist, and hence 
the construction mentioned above cannot be used to prove the integrability of the system. On the 
other hand, it turns out that not all (finite-dimensional) completely integrable systems admit a bi- 
Hamiltonian structure. In [B], R. Brouzet constructed a (four-dimensional) non-degenerate completely 
integrable system, (i.e., associated to a non-degenerate Poisson bracket,) which does not admit a bi- 
Hamiltonian structure. The geometry of non-degenerate completely integrable systems admitting 
a bi-Hamiltonian structure was later investigated in [BMT], and the results obtained in that paper 
strongly suggest that bi-Hamiltonian formulations for non-degenerate integrable systems will be rather 
exceptional. (To our knowledge, there is only one example of such a system, namely the Calogero 
system (see for example [MM]).) 

The problem of finding bi-Hamiltonian structures for degenerate completely integrable systems 
was investigated by many authors, and turns out to be completely different from the non-degenerate 
case. Indeed, a number of different techniques (e.g., i?-matrices, master symmetries, reductions of 
partial differential equations) were applied successfully to construct bi-Hamiltonian structures in this 
situation, and several important examples of (degenerate) completely integrable systems (e.g., the 
Toda lattice, some generalized Henon-Heiles potentials and the so-called master systems) were shown 
to admit bi-Hamiltonian structures. (We refer to [ARl],[AR2],[BR],[D],[Val] for examples.) However, 
a systematic procedure to find such a structure is, to our knowledge, not known at this moment. 

Starting from an arbitrary non-zero polynomial <p(x,y) G R[x,y], P. Vanhaecke [Va2] constructed 
a Poisson bracket {•, -} v on R M given, in terms of the standard coordinates {u\, . . . , Ud, vi,...,vo) by 



{u(X),u j r = {v(\),v j r = 0, 
{ulX;.c,}' K.riAlK ¥>(A,t;(A)) ^ " (A) 



\d+i-j 



mod u(X), 



for all j = 1, . . . , d, where u(x) = x d + + . . . + Ud and v(x) = vix d ~ 1 + . . . + Vd- (We have 

denoted by [i?(A)] + the polynomial part of the rational function R(X), and by P(A)mod u(X) the 
remainder of the division of P(X) by u(A).) Taking the product of this Poisson manifold with the 
trivial Poisson manifold (R ra+1 , {•, -} ) (with standard coordinates cq, . . . ,c n which, in what follows, 
will be identified with the coefficients of the polynomial c(x) = cqx 11 + . . . + c n ), one obtains a Poisson 
bracket {•, ■}<? on R M + n +! which has the functions Cq, . . . , c n as Casimir functions. 



Next, we consider deformation families of polynomials of the form 

F(x, y, c) = F(x, y) - x d c(x), c(x) = c x n + . . . + c n , 

where F(x, y) is a polynomial depending explicitly on y. It was shown in [Va2] that the (independent) 
coefficients Hi, . . . ,H d of the polynomial 

H(X) = HiX 1 ' 1 + ... + H d = (F(X, v(X)) - A d c(A))mod u(X) 

are in involution with respect to the Poisson bracket {•,■}£ (for every polynomial tp), and hence 
determine (together with the Casimir functions co,...,c n ,) a completely integrable system on the 
Poisson manifold (R 2d + n + 1 , {■, • }£). 

The first aim of this paper is to show that the functions in the set 

I = { c o, ■ ■ ■ ,c n , Hi, . . . , H d } 

actually Poisson-commute with respect to a much larger class of Poisson structures {•,■}£ (indexed by 
the polynomials <p and the integers a < n+1). Intuitively speaking, these new brackets are constructed 
from the Poisson bracket {•,■}£ by "choosing another set of Casimir functions". More precisely, the 
construction works as follows. For any integer a < n + 1, the polynomial 

G(X) = X d c(X) + {F{X,v(X)) - A d c(A))mod u{X) 

(which has the elements of / as its coefficients) can be written in the form 

G(X) = f X d+n + ... + f n - a X d+a + X a (g 1 X d - 1 + ...+g d ) + f n _ a+1 X a ~ l + ... + /„, 

(i.e., we make a special partition 

P = {{fo,---,fn},{9i,---,9d}} 

of the set I). The mapping 

<f> : Ti 2d+n+1 -> R M +»+! : ( n , Vj c ) _> ( u> W) f( u , v, c)) 

associated to the functions /o, ...,/„ is a bijection on (a dense subset M of) R M + n +! . Pushing forward 
the Poisson bracket {•, with the mapping we obtain a new Poisson bracket {•, (on M) which 
has the functions fo,---,f n as Casimir functions. The remaining (independent) functions gi, . . . ,g d 
are in involution with respect to the new Poisson bracket {•,•}£, and determine a (i-dimensional 
completely integrable system (M, {■, -}^,g)- 
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In the special case where <p = 1, a = and 

F(x,y)=y 2 + f(x), 

these completely integrable systems turn out to be the master systems introduced by D. Mumford 
[M] and P. Vanhaecke [Val], and all systems (M, {•, •}<£, g) constructed above will therefore be called 
generalized master systems. 

It is clear from the construction given above that (for a fixed polynomial F(x,y)) all generalized 
master systems are different in the sense that they are defined by different Poisson brackets and 
different Hamiltonians. However, it turns out that not all these systems are different in the sense 
that they have different Hamiltonian vector fields. To show this, we will compute the Hamiltonian 
vector fields of the generalized master systems, and we will search for different systems (i.e., different 
Poisson structures and Hamiltonians) yielding the same Hamiltonian vector fields. More precisely, we 
will consider the generalized master system (M, {■, -}^,<7), and we will search for other partitions P' 
(associated to the integer a' = a + i<n+l) and other polynomials <p' such that the Hamiltonian 
vector fields of the generalized master system (M, {•, ■}„/, 5') coincide with those of the first system, 
i.e., such that 

for all a = 1, . . . , d. We will show that this condition is satisfied when ip' = x l (p, thus proving that all 
Hamiltonian vector fields of the generalized master system (M, {■,■}%, g) can be written in Hamiltonian 
form with respect to the Poisson brackets {•, for all i such that a + i < n + 1. We will also 

show that all these Poisson brackets are compatible, leading to a multi-Hamiltonian formulation for 
"all" generalized master systems. (Strictly speaking, the above argument shows that our large class 
of generalized master systems consists of a (smaller) class of completely integrable systems (which 
turns out to be exactly the class introduced in [Va2] , if we allow 92 to be a special rational function) , 
together with the systems determined by their multi-Hamiltonian formulation, and from this (stricter) 
point of view, our construction yields a multi-Hamiltonian structure for the smaller class of systems 
introduced in [Va2].) 

Finally, we remark that the construction of the Poisson bracket {•,■}« can easily be adapted to 
the case where P is an arbitrary partition 

P = {{fo,---,fn},{9i,---,9d}} 

of the set /, yielding a new Poisson bracket which has the functions /o, ...,/« as Casimir 

functions. The remaining functions g 1 ,...,g c i again Poisson-commute with respect to this Poisson 
bracket and hence determine a completely integrable system, also called a generalized master system. 
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The procedure given above can be generalized to this larger class of completely integrable systems, 
again showing that our collection of systems consists of a (smaller) class of integrable systems, together 
with their multi-Hamiltonian formulation. (We remark however, that this class of (different) systems 
is larger than the one considered in [Va2], i.e., we have actually constructed new systems, in the sense 
that they have different Hamiltonian vector fields.) 

We have already remarked that, putting ip = 1, a = and F(x, y) = y 2 + f(x) in our description, 
we obtain the well-known (odd and even) master systems. In the two-dimensional case, the existence 
of a multi-Hamiltonian structure for these master systems was shown by P. Vanhaecke in [Val]. 
However, the construction given in [Val] involves a heuristic choice of the Casimir functions, and uses 
a trial-and-error method to compute the components of the different Poisson brackets involved in this 
multi-Hamiltonian structure. As an application of our results, we will give a complete description 
of the multi-Hamiltonian structure of the two-dimensional odd master system. We will compute the 
Hamiltonian vector fields of this system, and we will construct the different (compatible) Poisson 
brackets involved in its multi-Hamiltonian structure. 

The paper is organized as follows. In Section 2, we collect some preliminary material concerning 
bi-Hamiltonian structures and integrable systems, and we fix some notation which will be used in the 
rest of the paper. In Section 3 we start by briefly reviewing the results from [Va2]. Further, we will 
also introduce the generalized master systems, and we will prove the integrability of these systems. In 
Section 4, we describe the multi-Hamiltonian structure for these generalized master systems. Finally, 
in Section 5 we use the theory developed in this paper to construct a multi-Hamiltonian structure for 
the odd master system. 

2. Preliminaries 

A Poisson manifold (M, {•,•}) is a smooth manifold M, endowed with an anti-symmetric R- 
bilinear mapping 

{•, •} : C°°(M) x C°°(M) -» C°°(M) 

called a Poisson bracket, which is a derivation in both of its arguments, and which satisfies the Jacobi 
identity 

{{f,9},h} + {{g,h},f} + {{h,f},g} = 

for all f,g,h € C°°(M). Two functions f,g € C°°(M) are said to Poisson- commute, or to be in 
involution, if their Poisson bracket vanishes, i.e., 



{f,g} = o. 
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A function / which Poisson-commutes with all functions g G C°° (M) is called a Casimir function. As 
the Poisson bracket {•,■} is a derivation in its arguments, every function / G C°°(M) determines a 
vector field Xf on M, called the Hamiltonian vector field generated by /, by the relation 

X f (g) = {g,f} 

for all g G C°°(M). If / generates a non-zero vector field, (i.e., / is not a Casimir function,) / is said 
to be a Hamiltonian, and the triple (M, {•,•},/) is called a Hamiltonian system. 

Let (M, {•,•}) be a Poisson manifold of dimension n which admits k (independent) Casimir func- 
tions (i.e., it is of rank 2d = n — k). A Hamiltonian system (M, {■,•}, H) on this Poisson mani- 
fold is said to be Liouville integrable or completely integrable if (apart from the Casimir functions 
Ci,...,Cfc) it admits d independent invariant functions (sometimes called the Hamiltonians of the 
system) Hi = H,H 2 , ■ ■ ■ ,Hd, which are in involution with respect to the Poisson bracket {•, •}. We 
will frequently denote such an integrable system by (M, {•, •}, h), where h = (Hi, . . . , Hd) is a set of 
Hamiltonians of the system. 

Two Poisson brackets {•, -}i and {•, -} 2 on a manifold M are said to be compatible if their sum or, 
equivalently, any linear combination 

{-,-}i+A{-,-} 2 , 

is again a Poisson bracket. More generally, a collection of Poisson brackets are said to be compatible 
if any linear combination of them is again a Poisson bracket or, equivalently, if they are pairwise 
compatible. 

A Hamiltonian system (M, {•, -}i, Hi) is said to be multi- Hamiltonian if its Hamiltonian vector 

field 

= {■, Hi}i 

can be written in Hamiltonian form with respect to a set of compatible Poisson brackets, i.e., 

{■, Hi}i = {•, H 2 } 2 = ... = {•, H n } n 

for some functions H 2 , . . . , H n G C°°(M). In the special case where n = 2, we say that the system is 
bi- Hamiltonian. 

In the rest of this paper, we will use the following notational conventions. Let u(x) be a (generic) 
polynomial of degree d. Then (for a fixed integer a) any polynomial P(x) can be written, in a unique 
way, as 

P(x) = Q{x)u(x) + x a R{x), 

where deg R(x) < d. We will denote 

Q(x) = 
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P(x) 
u(x) 



Further, we define 

x a R(x) = P{x) - u{x) 

and we say that 

P{x) 
u(x) 

In the case where a = we will, in analogy with the notation in [Va2], denote 

P(x)modo u(x) = P(x)mod u(x), 

and 



~P(x)~ 




~P(x)~ 




~P(x)~ 




~P(x)~ 


u(x) 





u(x) 


+ 


u(x) 





u(x) 



3. Construction of generalized master systems 

In [Va2], P. Vanhaecke constructs, starting from any polynomial F(x,y) G R[^,y] \ R[z], a set 
of d independent functions Hi, . . . , H d . These functions turn out to be in involution with respect 
to a large class of Poisson brackets {•, on the space R M , and hence they determine a family of 
completely integrable Hamiltonian systems on this space. More generally, considering a deformation 
family 

F(x,y,c), c = (ci, . . . ,c fc ) G R fc , 
of polynomials, one can construct a new set 

I = {ci, . . . ,Cfc, Hi, . . . , 

of independent functions, which again turn out to Poisson-commute with respect to a family of Poisson 
brackets {•,•};? on R 2d + fc (which have the functions Ci,...,Ck as Casimir functions) and therefore 
determine a collection of completely integrable systems on the space R 2d+fc . The aim of this section is 
to show that the functions in the set I actually Poisson-commute with respect to a much larger class 
of Poisson brackets on (a dense subset M of) the space H 2d+k . Indeed, it will turn out that, choosing 
a "suitable" set fi, - ■ ■ ,fk of functions on the space R 2d+fc , we can construct a Poisson bracket {•, 
on a dense subset M C R 2d+fc , which has the functions fi,...,fk as Casimir functions. The functions 
in the set I are in involution with respect to this new bracket, and hence determine a (new) completely 
integrable system on M. 

We start this section by briefly reviewing the construction given in [Va2]. To this purpose, we 
consider the space R M as the space of pairs of polynomials 

U (A) = \ d +u 1 \ d ~ 1 + ... + u d , 

v{\)= viX*' 1 + . . . + v d , 
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P(x) 
u(x) 



P(x)mod a u(x), 



P(x) 




P(x) 


u(x) 


rv 


u(x) 



i.e., the coefficients of the polynomials u, v serve as coordinates on R 2d . The following result then 
gives explicit expressions (in terms of the coordinates iij, Vi on R 2d ) for a family of Poisson brackets, 
indexed by the space of (non-zero) polynomials in two variables: 

Proposition 1. Let (p(x,y) be a non-zero polynomial in R[x,y]. Then the bracket {-, given in 
terms of the coordinates Ui, Vi by 

{u(\), Uj r = {v(\), Vj y = o, 

{u(\),V j r = {u j ,vW} V =<P(Kv{\)) 



u(A) 



\d+l-j 



mod u(X), 



j = l,...,d, defines a (polynomial) Poisson structure on the space R 2d . Moreover, all these Poisson 
brackets {•, are compatible. 

The next step is the construction, for every Poisson bracket given by (1), of a collection of completely 
integrable systems on R 2d . This is the purpose of the following 

Proposition 2. Let F(x, y) G R[x, y] \ R[x]. Then the coefficients Hi, . . . , H& of the polynomial 

H(X) = H^- 1 + ... + H d = F(X, v(X))mod u(X) (2) 

define, for any non-zero polynomial (p(x,y), a completely integrable system on the Poisson manifold 
(R 2d , {•, with polynomial invariants. 



The previous construction can be extended to obtain completely integrable systems on certain 
product manifolds R 2d x D. We will now describe this procedure. However, we will restrict our 
attention to the case where D = R fc , as the more general situation will not be considered in the rest 
of this paper. 

We start by defining a (trivial) Poisson bracket {•, -}o on D = R fc by putting 

{/,<7}o = 

for all f,g G C 00 (R fe ). Letting (ci,...,Cfc) denote the standard coordinates on R fe and considering 
the space R 2d + fe as the product of the Poisson manifolds (R 2d , {-, and (R fc , {•, -}o), we obtain a 
(product) Poisson structure {•,•}£ on R 2d+fc , which has the coordinate functions c\, . . . , Ck as Casimir 
functions. Suppose that we are given a deformation family of polynomials, i.e., a family of polynomials 

F(x,y,c 1 ,...,c k ) 

depending on k parameters c = (ci, . . . , Ck) € R fc . (In the more general case, one takes parameters in 
the manifold D.) Then we have the following 
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Proposition 3. The coefficients H 1 , . . . ,Hd of the polynomial 

H(X) = F(X, v(X),ci, . . . , Cfc)mod u(X) 

are in involution with respect to the Poisson bracket {■, and hence determine a completely inte- 
grable Hamiltonian system (R 2d+k , {-, -}^, h), h = (H 1 , . . . , H d ) on the (degenerate) Poisson manifold 
(R M +M-,-}^). 

Let us suppose that (p(x,y) is a polynomial and that F(x,y,c), c = (ci, . . . , Cfc) £ R fc is a de- 
formation family of polynomials, and let (R 2d+fc , {•, ■}<?, h) be the completely integrable Hamiltonian 
system given by Proposition 3. We will now show how one can construct, starting from the system 
(R 2d+fe , {•, h), a large collection of other completely integrable systems on the space or, 
more precisely, on a dense subset M of H 2d+k ,) which have the same set of invariants. Roughly speak- 
ing, this construction works as follows: choosing a "suitable" set of functions /i, ...,/& in C 00 (R 2d+k ), 
we construct a new Poisson bracket {•,•}/ (on M) which has fi,---,fk as Casimir functions. We 
will show that the invariants ci, . . . ,Cfc, Hi, . . . , Hd of the original system (TL 2d+k , {•, ■}£, h) Poisson- 
commute with respect to {•, and therefore determine a completely integrable system. 

We start by constructing the new Poisson bracket {•, -}^. To this purpose, we consider k functions 

fi(c, h{u, v, c)), . . . , f k (c, h{u, v, c)) 

such that the mapping 

. Yi 2d+k -> R M+fc : ( u , «, c ) ' — ► (u, f(u, v, c)) 

(where we denoted / = (/i, . . . , /&)) is a bijection on a dense subset M C R M+fc . Pushing forward 
the (restriction to M of the) Poisson structure {•, •}<? with the mapping </> _1 , we obtain a new Pois- 
son bracket {•,■}/ on M. By construction, the mapping ^ is a Poisson morphism and hence the 
functions /j(c, h(u, v, c)), i = 1, . . . , k are Casimir functions for this new Poisson bracket. The follow- 
ing theorem states that the (restriction to M of the) invariants ci , . . . , c k , Hi , . . . , Hd of the system 
(R 2rf+fc , {■, •}£ , h) Poisson-commute with respect to this new bracket {•,■}?■ 

Theorem 4. Let fi(c, h(u, v, c)), . . . , fk(c, h(u, v, c)) be k functions such that the mapping 

4> : K 2d+k -> R M+fc : (u, u, c) ^ (u, u, /(it, v, c)) 

is bijective on a dense subset M C R 2d+fc , and Jet {•,■}/ be the Poisson bracket associated to f = 
(/i , . . . , /fc ) as above. Then the invariant functions Ci, . . . ,Ck, Hi, . . . , Hd of the system (R 2d+fc ,{-,-}^,h) 
are in involution with respect to the Poisson bracket {•,■}/ and hence determine a completely inte- 
grable Hamiltonian system on the Poisson manifold (M, {•, - }j). 
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Proof. For the sake of simplicity, we will denote u = (u\ , . . . , u^,, v\, . . . , Vd), and we write the Jacobian 
matrix J(f, c) of the functions fi with respect to the variables Cj by f c . 

As the functions ci, . . . , Ck are Casimir functions of the Poisson bracket {•,■}£, the Poisson matrix 
(with respect to the coordinates (u, c)) of the bracket {■, is of the form 

o oy ' 

where ^4 is the (skew-symmetric) (2d x 2d)-matrix given by {u, u} v . 

Next, we compute the Poisson matrix of the bracket {•, -}y with respect to the coordinates (u, c). 
As (f> is a Poisson morphism, 

{u,u}? = {u,u}£ = {u,u}*\ 
and hence its Poisson matrix is of the form 

A B\ 
- l B C ) ' 

Taking into account that the functions fi(c, h(u, c)), i = 1, . . . ,k are Casimir functions of the bracket 
{'>'}/! we obtain that 

J, *)(£)-a 

By assumption, det f c ^ on M and (3) yields 

B = -AU~\ C = t f- lt f u Af u f-\ (4) 

As the functions /j only depend on Cj and Hj (which Poisson-commute with respect to the bracket 
{•, ■}£), we see that 

{/,/}? = 7uA/„ = o, 

and hence 

C = 0. (5) 
Summarizing, the Poisson matrix of the bracket {-, is of the form 



^4 -^4/u/" 1 



(6) 



Using this expression, we compute that 
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and as all the functions fa and Hj are in involution with respect to {■,•}£, this proves that the 
Hi Poisson-commute with respect to Analogous computations yield that {c,h}^ = and 

{c, c}j = 0. Hence, we obtain (apart from the Casimir functions fi,---,fk) a set of d additional 
independent functions in involution, which concludes the proof of the theorem. 

I 

In the rest of this section, we discuss an application of Theorem 4 which will play an important 
role in what follows. We will consider special deformation families of polynomials. Inspired by our 
study of the master systems (described in Section 5), we make a very simple and natural choice for 
the functions / = • • • , /jt). It turns out that the mapping <p associated to / is bijective on a 
dense subset M C R 2d+fc , and applying Theorem 4 we obtain a completely integrable system, called 
a generalized master system, on the Poisson manifold (M, {•, ■}%)■ 

We start by fixing an integer d > 1, a (non-zero) polynomial ip(x,y) and a polynomial 

F(x,y) €K[x,y]\R[x], 

and we consider the deformation family 

F(x,y,c ,...,c n ) = F{x,y) -x d c{x), (c , . . . , c n ) G R n+1 , (7) 

where c(x) = cqx 71 + . . . + c n . Further, we denote by (R 2d+ra+1 , {-, •}£, h) the completely integrable 
system constructed in Proposition 3, and by I the set 

I = { c 0, ■ ■ ■ ,Cn, Hi, . . . , H d } 

of invariants of this system. Finally, we also introduce the polynomial 

G(X) = X d c(X) + H(X) = X d c(X) + (F(X, v(X)) - A d c(A))mod u{X), (8) 

whose coefficients are given by the elements of the set /. 

Choosing any integer a < n + 1, we can write the polynomial G(X) in the form 

G(X) = f X d+n + ... + f n - a X d+a + X a (g l X d ~ 1 + ...+g d ) + U.^X*- 1 + ... + /„, (9) 

and we obtain a partition 

P = {{fo,---,fn},{9i,---,9d}} 

of the set / = {co, . . . , c n , Hi, . . . , Hd}. Our aim is to apply Theorem 4, where / = (fo, ■ ■ ■ , f n ) is 
given by the partition P (or, equivalently, by the integer a). It is clear that, in order to do this, we 
have to show that the mapping 

. R 2d+n+l ^ R 2d+n+l . ^ ^ ^ ^ ^ y ^ y ^ ^ 
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associated to / = (f , . . . , f n ) is bijective on a dense subset M C R 2d + ra + 1 . 
It is easy to see that (8) implies 

A d c(A)" 



u(X) 



G(A) = F(X, v(X))mod u(X) + u(X) 
On the other hand, denoting 

fW = foX d+n + . . . + f n - a X d+a + f n - a+ iX a 1 + . . . f n , 

g(X) =g 1 X d ~ 1 + . . . + g d , 

(9) takes the form 

G(X) = f(X) + X a g(X). 

Comparing (10) and (11) we obtain 



which yields 



(10) 



(11) 



"A rf c(A)~ 




r/(A)-(F(A,t;(A))mod«(A))- 






+ 




a 



c(A) 



u(A) 
A d 



/(A)-(F(A, V (A))mod u(X)) 



u(X) 



(12) 



The right hand side of this equation is defined for a generic polynomial u(A). More precisely, it is 
defined for all polynomials u(X) in the case where a = 0, and on the subspace u& ^ in the other 
cases. Hence, the mapping 4> is bijective on the dense subset M = H 2d+n+1 (\{ud = 0}). As an 
immediate consequence of Theorem 4, we obtain the following 

Corollary 5. Let d > 1 be an integer and ip(x,y) any non-zero polynomial in H[x,y]. Further, let 

F(x, y, c) = F(x, y) - x d c(x), c(x) = c x n + . . . + c n , 
be a deformation family of polynomials and denote by 

G(X) = X d c(X) + (F(X,v(X)) - A d c(A))mod u{X) 
the polynomial whose coefficients are given by the set 

I = { c 0, ■ ■ ■ ,Cn, Hi, . . . , H d } 

of invariants of the system (H 2d+n+1 , {•, h). For every a < n + 1, the partition 

P = {{fo, ■ ■ ■ , fn},{9i, ■ ■ ■ ,9d}} 

of the set I, given by 

G{X) = f X d+n + ... + f n - a X d+a + X a (g l X d ~ 1 + ... + g d ) + f n _ a+1 X a - 1 + ... + /„, 

determines a completely integrable Hamiltonian system (M,{-,-}^,g) on the dense subset M = 
p2d+n+i^j n ^ _ called a generalized master system. 
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4. Multi-Hamiltonian formulation for the generalized master systems 

The main purpose of this section is to construct a multi-Hamiltonian formulation for the gen- 
eralized master systems (M, {■,•}„,<?) introduced in Corollary 5. More precisely, we will prove (the 
slightly stronger result) that all Hamiltonian vector fields X ga , a = l,...,d, (not only X gi ) of the 
generalized master system (M, {•, •}<£, g) can be written in Hamiltonian form with respect to a set of 
compatible Poisson brackets. The first step in our construction is to determine explicit expressions 
for the Hamiltonian vector fields X 9a = {■, <7 a }a, a = l,...,d, of the generalized master system with 
respect to the coordinates (u,v,c). Next, we show that the Hamiltonian vector fields 

x ga={ m >9a}%, a=l,...,d, 
of this system coincide with the Hamiltonian vector fields of a number of other (generalized master) 
systems (M, {•, g') (for particular choices of the polynomial ip' and of the integer a'), i.e., 

{') 9a} a = {') 9 a }a' ' 

for all a = 1, . . . , d. As an immediate consequence of this fact, the Hamiltonian vector fields of the 
generalized master system (M, {-,■}«> g) can be written in Hamiltonian form with respect to a set of 
Poisson brackets {-, -}^, (for certain (p' and a'). Finally, we will show that any pair of these Poisson 
brackets are compatible, thus concluding our construction of a multi-Hamiltonian formulation for the 
generalized master systems. 

We start our construction by computing explicit expressions for the Hamiltonian vector fields 
of the generalized master system (M, {■,■}%, g) with respect to the coordinates (u, v,c). This is the 
purpose of the following 

Theorem 6. Let F(x, y, c) be a deformation family of polynomials as above. Further, suppose that 
a < n + 1 is an integer and let 

P = {{fo,---,fn},{9i,---,9d}} 

be the partition of the set 

I = { c o, ■ ■ ■ ,c n , Hi, . . . , H d } 
associated to a, i.e., the partition given by 

G(X) = f X d+n + ... + f n - a \ d+a + A^iA^" 1 + . . . + g d ) + /n-^iA"" 1 + ... + /„. (13) 

The Hamiltonian vector fields X ga = {-,g a }a °f tne generalized master system (M, {■, g) are 

completely determined by (the coefficients in /j, of) 

dF 

X G (v)u(\) = -^j{fi,v(fi)){u(X),v(n)} v mod a u(fi), 

-F(,M,))-^c W MX)Mfl)VmodaU{fiy (14) 



u(n) 



+ 
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Proof. As the functions fo,---,f n are Casimir functions of the Poisson bracket {•, •}£, (13) yields 

{, G(»)Z = ^({-^ilS/- 1 + • • • + {;9d)t)- (15) 



On the other hand, we can rewrite G(/j,) as 
which implies that 

,dF 



u(p) 



J + 



Comparing (15) and (16), we obtain 



F(p,v(^) - ^ d c(n) 
u(n) 



^(^-Lx gi + ...+x gd ) = ({-,v(p,)r—^,v(ri) 



BF 



dy 



u(n) 



mod Q u((x), 



which proves (14). 

Finally, substituting 



in the expression 



G(A) 



A d c(A) 



u(X)+F(X,v(X))mod u(X) 



{G(X),g a Z = 0, 



(16) 



we see that, for all a = 1, . . . , d, X ga c{X) is completely determined by X ga u(X) and X ga v(X), which 
concludes the proof of the theorem. 



Our next step is to search for polynomials ip' and partitions P' (associated to a'), such that the 
Hamiltonian vector fields 

x ga = {-,9a}%, a = l,...,d, 

of the generalized master system (M, {•, g) coincide with the Hamiltonian vector fields of the 
system (M, {•, ■}„,,</'), i.e., 

for all a = 1, . . . , d. 

We start by stating a simple lemma, which will turn out to be very useful. 
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Lemma 7. Let a and i be integers, P{x) an arbitrary polynomial, and u(x) = x d + u\x d 1 + . . . + Ud 
a generic polynomial. Then 



x l P(x) 



u\x) 



P x) 



u\x) 



and 



(x l P(x))mod Q+ i u(x) = x t (P(x)mod a u(x)). 



(17) 
(18) 



Proof. If 

P(x) =u(x)Q(x) + x a R(x), 

where degR(x) < deg-u(x), then 

x l P{x) = u{x)x i Q{x) + x a+i R(x), 
and the result follows immediately. 

I 

Let ip, F(x, y, c) and G(X) be as above and let a < n + 1 be an integer. Further, we denote by 

P = {{fo,---,fn},{9i,---,9d}} 

the partition of the set / associated to a and we consider the generalized master system (M, {•, -}%,g) 
associated to this partition. It is clear from Theorem 6 that 



H a (fi d ~ 1 X gi u(\) + ... + X gd u{^)) = -^{n.,v{n)){u{\),v{v)yraod a u(fi). (19) 



Let iff be another polynomial, suppose i is an integer such that a + i < n + 1 and denote by 

the partition of the set /, associated to a' = a + i. Applying Theorem 6 for the generalized master 
system (M, {-, -}^ +i ,g'), we obtain 

/j, a+i (fi d ~ 1 X' g ^u(\) + ... + X' g ,u{\)) = —{jj., «(m)){«(A), v(tx)r'mod a+i u(»). (20) 

Our aim is to search for a polynomial tp' such that the Hamiltonian vector fields X 9a and X' gl coincide 
for all a = 1, . . . , d. Comparing (19) and (20), we obtain the (necessary) condition 

/z l ({u(A), w(/x)}^mod a u(n)) = {u(A), v (fi)}^ mod a+i u(n), (21) 
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which, in view of Lemma 7, is equivalent to 

^{u(X),v(fi)} (p mod a+i u(n) = MA),u(/x)} v mod a+i u(fi). 
Using Proposition 1, we see that, for a = 1, . . . , d, 

a 1 {u a ,v{ij)Y m.od a+i u{n) =n' l {u(fi),v a } ip mod a+ i u{n) 



(22) 



,,d+l — a 

u(n) 



d+l-a 



=(n l <p{n,v(n)) 

={u(p),v a } xlip mod a+ i u(/x) 
={u a ,v(fi)} xtip mod a+ i u(n), 



mod it(/u) )mod a+i u(fi) 



IP 



d+l-a 



mod a+i u(ii) 
mod u(n))mod a+ i u(/i) 



J + 



(23) 



and hence, 

fi l {u(X),v(p)} ip mod a+% u(p) = {u(X), v (//)} xV mod a+i u(/x). (24) 

Comparison of the expressions (22) and (24) suggests taking ip' = x % ip, and we obtain the following 
Theorem 8. Let a and i be integers such that a + i < n + 1, and let 

P = {{f ,...,f n },{ gi ,...,g d }} and P' = {{/£, . . . , f n }, {g[, . . . , g' d }} 

be the partitions of the set I, associated to a and a' = a + i. Then 

^({«(A), v(u,)ymod a u(jj)) = MA), v{n)Y^mod a+i u(ji), (25) 

and the Hamiltonian vector fields of the generalized master systems (M, {■ , •}£ , g) and (M, { ■ , • j d' ) 
coincide, i.e., 

for all a = 1, . . . , d. 



Proof. The first statement follows immediately from (24) and Lemma 7. Using Theorem 6 and (25), 
we see that 

{«(A),<£}fe={«(A),aJS, , , 

(26) 

MA),fli}&=MA),ft,}S, 

for all a = 1, . . . , d. As the Hamiltonian vector fields X 9a and X' gl are completely determined by (14), 
this proves the second part of the theorem. 
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As a consequence of this theorem, we can write the Hamiltonian vector fields of the generalized 
master system (M, {•, •}<£, g) in Hamiltonian form with respect to a set of Poisson brackets {•, -}^_^, 
indexed by the integers i such that a + i < n + 1. We will now show that all these Poisson brackets 
are compatible. To this purpose, it suffices to show that the Poisson brackets are pairwise compatible, 
i.e., that for any polynomial 92 and any integers a and i such that a + i < n + 1, the Poisson brackets 
{■,•}<£ and {•, -}q_J^ are compatible. To prove this assertion, we will use a very interesting trick which 
was, to our knowledge, first used by P. Vanhaecke in [Val], to show that the Poisson brackets of a two- 
dimensional completely bi-Hamiltonian integrable system (we refer to [Val] for the exact definition of 
such systems) are compatible. 

Theorem 9. Let ip be an arbitrary polynomial, and suppose that a and i are integers such that 
a + i < n + 1 . Further, we denote by 

P = {{fo,...J n },{9u---,9d}} and P' = {{/£, . . . , f' n }, {g[, . . . , g' d }} 

the partitions of the set I associated to a and a' = a + i. Then the Poisson brackets {•, ■}£ and 
{■, '}a+i °f the generalized master systems (M, {•, •}<£, g) and (M, {•, ■} x ol ^ i -,g') are compatible. 

Proof. To simplify the notation, we denote {•, -}i = {•, ■}£ and {-, -} 2 = {•, Further, we write 

the polynomial G(X) as 

G(A) = H X d+n + ... + H d+n , 

and we denote by 

{v} = {v}i + {v} a 

the sum of the two brackets. 

To prove the compatibility of {•, -}i and {•, -}2, it suffices to prove the Jacobi identity 

{{/, g}, h} + {{g, h}, /} + {{h, /}, g} = (27) 

for all functions f,g,h G {ui, . . . , Ud, vi,...,Vd, H 0l . . . , Hd+n}- Using the fact that {•, -}i and {•, 
are Poisson brackets and hence satisfy the Jacobi identity, (27) reduces to the "mixed Jacobi identity" 

{{/, gh,hh + {{g, hhJh + {{h, fh,gh 

(28) 

+ {{/, 9h , hh + {{9, hh , f} 2 + {{h, />! , g} 2 = 0. 
From Proposition 1 we know that {•, -} v and {•, -} xtip are compatible. Hence, (28) holds if f,g,h € 
{ui, . . . , Ud, v\, . . . , Vd}, and we only have to prove (28) in the case where at least one function (say, 
h,) belongs to {H , . . . , Hd+ n }- As Theorem 8 implies that, for alH = 0, . . . , d + n, we have 

{;H i } 2 = {;H j } 1 (29) 
{;Hi}i = {;H k } 2 (30) 
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for some j and k, we see that 

{{f,g} 2 , Hi}! + {{g, Hi} 2 , f}, + {{Hi, f} 2 ,g}i 
+ {{f,g}i,Hi} 2 + {{g,Hi} 1 ,f} 2 + {{Hi,f} 1 ,g} 2 
={{f,g} 2 ,H k } 2 + {{g, Hj} u fh + {{Hj, fh^h 

+ {{/, gyuHjh + {{g, H k } 2 ,f} 2 + {{H k , f} 2 , g} 2 

=0, 

as {•, -}i and {•, -} 2 satisfy the Jacobi identity. 

I 

Combining the results stated in this section, we obtain the following 

Theorem 10. Let F(x, y, c) = F(x, y) — x d c(x) be a deformation family of polynomials, 92 an 
arbitrary non-zero polynomial, and 

G{\) = \ d c{\) + (F(A,u(A)) - A d c(A))mod u{\) 

the polynomial whose coefficients are given by the set 

I = { c o, ■ ■ ■ ,c n , Hi, . . . , H d } 

of invariants of the system (R M+n+1 , {•, ■}£, h). Further, suppose that a < n + 1 is an integer and let 

P = {{fo,---,fn},{9i,---,9d}} 

be the partition of the set I associated to a, i.e., the partition given by 

G(X) = f \ d+n + ... + f n - a X d+a + A"^- 1 + ...+g d ) + f n - a+ i\ a - 1 + ... + /„. 

The generalized master system (M, {■,•}<£, 3) admits a multi-Hamiltonian structure, which is con- 
structed as follows. For every integer a' = a + i<n + l, the partition 

P' = {{f[ ) ,...,a,{g[,...,g' d }} 

of the set I, associated to a', determines a Poisson bracket {•, -}^_^, and the Hamiltonian vector fields 
{-,g a }a> a = 1) ■ ■ ■ 5 rf> of the system can be written in Hamiltonian form 

with respect to the (compatible) Poisson bracket {■, - j^+i- 
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Remark 1. Strictly speaking, it follows from Theorem 8 that well-chosen families of generalized 
master systems have the same Hamiltonian vector fields (and will therefore not be considered as 
different systems in the strict sense). Combining this result with Theorem 9, we see that the class 
of generalized master systems therefore consists of a (smaller) class of completely integrable systems, 
together with the systems determined by their multi-Hamiltonian structure, and from this stricter 
point of view, the construction given in our paper actually yields a multi-Hamiltonian formulation 
for a smaller class of completely integrable systems (instead of, as was stated in Theorem 10, for all 
generalized master systems). 

Remark 2. In Corollary 5, we have constructed a class of generalized master systems, starting from 
a (special) partition 

P = {{fo,---,fn},{9i,---,9d}} 
of the set /, associated to an integer a < n + 1, i.e., determined by 

G(X) = f X d+n + ... + f n - a \ d+a + A^iA^" 1 + . . . + g d ) + f n - a+ i\ a ~ l + ... + /„. 
Choosing an arbitrary partition 

P = {{fo, ■ ■ ■ , fn},{9l, ■ ■ ■ ,9d}} 

of the set /, the mapping (f> (associated to /) is again bijective on a dense set M in R 2d+fc , and 
applying Theorem 4 in this case, we again obtain a completely integrable system on M, which will 
also be called a generalized master system. In the rest of this section, we will show how to adapt our 
construction of a multi-Hamiltonian formulation to the case of arbitrary partitions. 

We start by introducing some terminology, similar to that in Section 2. Let P(x) be an arbitrary 
polynomial and assume that u(x) = x d + uix d_1 + . . . + u& is a generic polynomial. Further, let 
a = (ai, . . . , ad) be a fixed d-tuple of different integers. The polynomial P(x) can be written, in a 
unique way, as 

d 

P(x) = Q(x)u(x) + ^R iX a \ 

i=l 

and as in Section 2 we use the notation 

d 

R(x) = ^^RiX ai = P(x)raod a u(x). 

i=i 

An arbitrary partition 

P = {{fo,---,fn},{9i,---,9d}} 

of the set 

I = {co, ■ ■ ■ ,Cn, Hi, . . . , H d } 
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Q(x) 



P{x) 
u(x) 



corresponds to a ci-tuple a = (ai, . . . , ad), gj being the coefficient of x a i in G(X). (We remark that this 
d-tuple satisfies the additional condition that a < d + n, i.e., a,j < d + n for all j = 1, . . . , d.) Using 
this correspondence between partitions and d-tuples, together with the notation introduced before, 
the expressions for the Hamiltonian vector fields given in Theorem 6 remain valid in the case of a 
general partition. 

To construct the multi-Hamiltonian formulation in this more general case, we consider two par- 
titions 

P = {{fo,...,fn},{9i,---,9d}} and P> = {{f^...,ti},{g[,...,g d }} 

of the set I, corresponding to two ti-tuples a = (ai, . . . , ad) and a + i = (ai + i, . . . , ad + i)- One 
can then show (by the same argument as before) that the Hamiltonian vector fields of the generalized 
master systems (M, {•, ,g) and (M, {•, •} X f,' p ,g') coincide, and that the Poisson brackets {•, -}j and 
{•, - } X fi V are again compatible, which concludes our construction. 

5. Multi-Hamiltonian formulation for the odd master system 

This last section is devoted to the study of a simple but very illustrative example of a (two- 
dimensional) generalized master system. More precisely, we consider the completely integrable system 
(R 7 , {•, -}(),g), associated to the polynomial 

F(x,y) = y 2 + f(x), deg/(x) = 5. 

We compute the Hamiltonian vector fields of this system, and we construct a multi-Hamiltonian 
formulation for this system by using the theory developed in the previous section. 

We start our investigation by making some notational conventions. We put d = 2, n = 2, and we 
denote 

u(X) = X 2 +u 1 X + u 2 , 

v(X) = uiA + v 2 , 

c(A) = c X 2 +c 1 X + c 2 . 
Further, we fix a polynomial F(x,y), given by 

F(x, y) = f{x) + y 2 , f{x) = x 5 + Ax i + Bx 3 + Cx 2 + Dx + E. 

Finally, for reasons that will become clear at the end of this section, we replace the variables c = 
(co,ci,C2) by a new set of variables (w\, W2, ws), given by 

w(X) =A 3 + w\X 2 + W2X + W3 
_ ~ F(X,v{X)) - X 2 c{X) ' 

~ _ ^(x) J + - 
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In these new variables, the polynomial 



G(X) =A 2 c(A) + {F(X,v{X)) - A 2 c(A))mod u{\) 
=F(X,v(X)) -u(X)w(X) 



takes the form 



G(X) = H^ 4 + H 2 X 3 + H 3 X 2 + H 4 X + H 5 , 



where 



Hi = A — m — w\ , 

H 2 = B - u 2 - u\W\ - w 2 , 

H 3 = C + v\ - u 2 wi - u 1 w 2 - w 3 , 

H 4 = D + 2vtv 2 - u 2 w 2 - utws, 

H 5 = E + v\ - u 2 w 3 . 

We will now compute the Poisson bracket {•, -}J associated to the partition 

Pn = \ f = 



>o = {/ = {H u H 2 , H 3 }, g = {H 4 , H 5 }}. 



Putting ip = 1 in Proposition 1, we see that the Poisson matrix of the bracket {•, -} 1 (on R 4 ) is of the 
form 

/ 1 \ 

A = 



1 ui 
-10 
V-l -ui o/ 



To compute the remaining components of the Poisson bracket {-,-}o, we use the following simple 
method: requiring Hi, H 2 and H 3 to be Casimir functions of the Poisson bracket, i.e., {-,Hi}l = 
for i = 1,2,3, we obtain a system of equations in the components of the Poisson bracket, which is 
then solved for the unknown components. Using this method, we find that the Poisson bracket {•, -}q 
is given (with respect to the coordinates (u,v,w)) by 



( 











1 

















1 


Ul 








2vi 







-1 











1 


Wl — Ui 




-1 


-Ui 








1 


Wl 


w 2 - u 2 













-1 











\ 








-1 


-Wl 








-2vi 





-2vi 


Ui — Wi 


u 2 - w 2 





2vi 






(31) 



The Hamiltonian vector fields Xh 4 and Xh 5 of the system (R 7 , {•, - }l,g) take the following form 
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(this can be computed directly from (31), or by using Theorem 6): 



' Ul 


= 2vi, 






= 2v 2 , 


U 2 


= 2t> 2 , 






= 2v 2 u 1 - 2v 1 u 2 , 


Vl 


= W 2 - U 2 - UtW! + u\, 






= W 3 - U 2 Wi + U\U2, 


< V 2 


= W 3 - W 2 Ui + U\U2, 






= uiw 3 - w 2 u 2 + ul, 


W\ 


= -2«i, 






= ~2v 2 , 


w 2 


= — 2w 2 — 2t>itfi + 2v\Ui, 




w 2 


= -2v 2 wi + 2v\u 2 , 


. ^3 


= 2uiv 2 - 2v 2 w 1 , 






= 2viwz + 2u 2 v 2 - 2v 2 w 2 



Using the theory developed in Section 4, the Hamiltonian vector fields (32) of this system can be 
written in Hamiltonian form with respect to three additional (compatible) Poisson brackets {•, -}f\ 
i = 1, 2, 3. We will now compute these brackets explicitly. 

First, putting ip = x in Proposition 1, the Poisson matrix of the bracket {•, -} x takes the form 

/ 1 \ 
-u 2 I 

-100 or 

V u 2 / 

Using the fact that Hi, H 2 and H§ are Casimir functions of the bracket its Poisson matrix 

takes the form 



( 





1 






















-u 2 








-2v 2 


-1 











1 


Wl — Ui 








u 2 











-u 2 


-w 3 








-1 




















Ui — Wi 


u 2 








2v 2 


V 


2v 2 





w 3 





-2v 2 


/ 



and it is clear (from Theorem 8 or by direct computation) that Xh 4 = {-,H s }f and Xh 5 = { m jH^}f, 
showing the Hamiltonian form of the Hamiltonian vector fields (32) with respect to this second bracket. 

2 

A similar computation yields the following Poisson matrix for the Poisson bracket {•, -} 2 : 



( 





-Ui 


-ui 





-2vi 


-2v 2 


\ 








-u 2 








-2v 2 





Ui 


u 2 








-Ui 


-w 2 


-U>3 




u 2 











-u 2 


-w 3 













Ui 


u 2 





2vi 


2v 2 




2vi 


2v 2 


w 2 


w 3 


-2vi 








/ 


\2v 2 





w 3 





-2v 2 









(Again, it is easily checked that Xh 4 and can be written in Hamiltonian form with respect to 
this new bracket.) Finally, putting ip = x 3 in Proposition 1, we see that {•, -} x3 is given by the matrix 

/ uj - u 2 uiu 2 \ 

A= uiu 2 uj I 

—uf + u 2 —uiu 2 I 

V -uiu 2 -ul / 
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Requiring H 3 , H4 and H 5 to be the Casimir functions of the bracket {•, -}| , we see that its Poisson 
matrix is of the form 

( A B ) 

where 

/ — 2v\ 2u\V\ — 2v 2 2v 2 u 1 \ 
■q _ -2v 2 2u 2 v 1 2v 2 u 2 

-W 2 + W\U\ ~W 3 + U 1 W 2 U1W3 ' 
V — 10 3 + W\U 2 u 2 w 2 u 2 w 3 I 

and 

(0 — 2v\W\ —2v 2 w\ \ 

2t;iu;i 2wittJ 3 — 2v 2 w 2 , 

2v 2 w\ 2v 2 w 2 — 2viw 3 y 

and we can again write the vector fields (32) in Hamiltonian form with respect to this bracket. 

Remark 3. In [M], D. Mumford constructs a natural vector field on the Jacobian of any hyperelliptic 
curve of genus g, given by an equation y 2 = f(x), where deg f(x) = 2g + l. This vector field turns out 
to determine a g-dimensional completely integrable Hamiltonian system. The construction of these 
integrable systems was later adapted to the case where deg/(x) = 2g + 2 by P. Vanhaecke [Val], who 
also made a detailed study of these systems, which he called the odd (resp. even) master systems, in 
the two-dimensional case (i.e., g = 2). 

The existence of a multi-Hamiltonian formulation for the two-dimensional master systems was 
already shown by P. Vanhaecke [Val]. However, his construction was of a heuristic nature, and the 
computation of the compatible Poisson brackets required a trial-and-error method. 

Comparing the differential equations (32) to the expressions given in [Val], we see that the two- 
dimensional odd master system and the system (R 7 , {•, -}q, g) given by (32) are the same (up to a 
constant factor), which provides us with a (more systematic) construction of the multi-Hamiltonian 
formulation for the (two-dimensional) odd master system given in [Val]. A similar computation, 
starting from the polynomial 

F(x, y) = y 2 + x 6 + Ax 5 + Bx 4 + Cx 3 + Dx 2 + Ex + F, 

yields a multi-Hamiltonian formulation for the (two-dimensional) even master system. 
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